format short e
echo on
% definizione della funzione 
F = inline('1-(1-x).^(0.5)'); 
pause
% estremi dellíntervallo di integrazione
a=0;
b=1; 
pause
% valore esatto dell' integrale
Iteor = 1/3;
pause
% fomula dei trapezi
It = trapezi(F,a,b)
pause
% formula di Simpson
Is = simpson(F,a,b)
pause
% formula dei trapezi composta
Itc = ctrap(F,a,b,11)
pause
% aumentiamo n
Itc = ctrap(F,a,b,21)
pause
% calcolo dell'integrale con errore minore di tol
tol = 1e-4;
pause
% calcolo dell'integrale usando la formula di Simposon adattativa
[Isa, fis] = quad(F,a,b,tol)
pause
% traccia dell'algoritmo [fnc,a,b-a,Is]
[Isa, fis] = quad(F,a,b,tol,1);
 
